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Abstract 

The analytical nonautonomous rogons are reported for the inhomogeneous nonlinear Schrodinger equation 
with variable coefficients in terms of rational-like functions by using the similarity transformation and direct 
ansatz. These obtained solutions can be used to describe the possible formation mechanisms for optical, 
oceanic, and matter rogue wave phenomenon in optical fibres, the deep ocean, and Bose-Einstein condensates, 
respectively. Moreover, the snake propagation traces and the fascinating interactions of two nonautonomous 
rogons are generated for the chosen different parameters. The obtained nonautonomous rogons may excite 
the possibility of relative experiments and potential applications for the rogue wave phenomenon in the field 
of nonlinear science. 
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The nonlinear Schrodinger (NLS) equation is a foundational model describing numerous nonlinear physical 
phenomenon in the field of nonlinear science such as optical solitons in optical fibres [TJ [21 [3] , solitons in the 
mean- field theory of Bose-Einstein condensates [4j [5l E] , and the rogue waves (RWs) (also known as freak waves, 
monster waves, killer waves, giant waves or extreme waves) in the nonlinear oceanography [7J [El [HI [TO] • RWs 
are single waves generated in the ocean with amplitudes much higher than the average wave crests around 
them [TT] [T3]. Recently, the RWs have attracted more and more attention from the point views of both 
theoretical analysis [21j[2l[23l[2l[25l[28l[29l[34]and experimental realization [13l HH [30l [31j [32l [33] . The oceanic 
RWs can be, under the nonlinear theories of ocean waves, modelled by the dimensionless NLS equation [8] 



.dtp ld 2 tp 



i- 



|VtV = o, (i) 



dt 2dx 2 

which describes the two-dimensional quasi-periodic deep-water trains in the lowest order in wave steepness and 
spectral width. In addition, it has been shown that the RWs can be generated in nonlinear optical systems 
and the term "optical rogue waves" was coined by observing optical pulse propagation in the generalized NLS 
equaiton [T3JIII]- Some types of exact solutions of Eq. ([T]) have been presented to describe the possible formation 
mechanisms for the RW phenomenon such as the algebraic breathers (Peregrine solitons) [15] , the time periodic 
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breather (Ma solitons) [16] , the space periodic breathers (Akhmediev breathers) \F7\ [18] . More recently, the 
Akhmediev breathers were also further studied [191 I2Ql EU [22] ■ A possible mechanism for the formation of RWs 
was also exhibited by using two-dimensional coupled NLS equations [241 1251 126) describing the nonlinearly 
interacting two-dimensional waves in deep water. The three-dimensional mechanism of RW formation was 
studied in a late stage of the modulational instability of a perturbed Stokes deep-water wave [27] . Furthermore, 
the optical RWs were also found in the NLS equation with perturbing terms (the higher-order NLS equation 
with the third-order dispersion, self-steeping, and self- frequency shift) |23) . Recently, the existence of matter 
RWs in Bose-Einstein condensates was predicted to cither load into a parabolic trap or embed in an optical 
lattice |34) . 

Here we point out that based on physical similarities between the rogue waves and solitary waves, hrst 
observed by Russell in 1834 and further known as solitons by Zabusky and Kruskal in 1965 [45] . we coin the 
word "Rogon" for each such "Rogue Wave" (or the word "Freakon" for the "Freak Wave") if they reappear 
virtually unaffected in size or shape shortly after their interactions. Similarly, there also exist the corresponding 
new terms "oceanic rogons" , "optical rogons" , and "matter rogons" for the oceanic rogue waves [71 [H [9[ [10] , 
optical rogue waves [T3l [14l [2T] [22j [23 , 28, 29], and matter rogue waves [34], respectively. 

To the best of our knowledge, there was no report on exact solutions related to the modified RWs (rogons) 
with variable functions before. In this Letter, we will extend the NLS equation ([1} to the inhomogeneous NLS 
equation with variable coefficients, including group- velocity dispersion /3(f), linear potential v(x,t), nonlinearity 
g(t) and the gain/loss term j(t), in the form [35[ [36l [37j [38] 

*^ + + W + ffWMV = nm (2) 

which is associated with 8£/5ip* — 0, in which the Lagrangian density is written as C = i(il>tip* — V'V't) — 
f3(t)\tfj x \ 2 + g(t)\ip\ 4 + 2[v(x,t) — i-f(t)]\ip\ 2 , where tp = ip(x,t), and tp* denotes the complex conjugate of the 
physical field ip. Eq. @ can be also known as the generalized Gross-Pitaevskii equation with variable coefficients 
for /3(t) = 1 [U [5[ [S[ [33]. Our goal is focused on the first-order and second-order rational-like solutions of Eq. 
@ to describe the possible formation mechanisms of optical RWs by using the similarity transformations [39, 
|40l [3U [42j [43j |44] and direct ansatz [21] [22]. Moreover, we analyze the dynamical behaviors of these solutions 
and interactions of two optical RWs (rogons) by choosing the different functions. 

Similarity transformation and nonautonomous rogons — To investigate the analytical rational-like solutions 
of Eq.@ related to the optical nonautonomous rogons, we employ the envelope field ip(x,t) in the gauge 

form gU [321 S3 S3] 

ip(x, t) = [* R (x, t) + i*i(x, t)} eW*'*) (3) 

whose intensity can be written as \ip(x,t)\ 2 = \^ji(x, t)\ 2 + \*f?i(x, t)\ 2 , where ^B,{x,t), tyj(x,t) and <p(x,t) are 
real functions of space-time x, t. The substitution of transformation ([3]) into Eq. ([2]) yields the system of coupled 
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real partial differential equations with variable coefficients 

*R,t + ^ (*i,xx + 2ip x ^ RtX - tpl^i + tp xx V R ) + [v(x, t) - ip t ]Vi + g(t)(^l + - 7(t)* R - 0, (4a) 

-*i,t + (*r, m - 2^ x *i >a! - y£*R " + *) " Vt]* r + + *l)*R + 7(*)*i = °- ( 4b ) 

By introducing the new variables rj(x, t) and 7~(f), we further utilize the following similarity transformations 
for the real functions ^> R (x,t), ^j(x,t) and the phase <p{x,t) 

* R (x, t) = A(i) + B(t)P( V (x, t), r{t)), 

9i{x,t)=C(t)Q(ri(x,t) > T(t)), (5) 
<Pi%,t) = x(x,t) + fir(t), 

to system (j4]) such that we deduce the following similarity reduction 

Vxx = 0, (6a) 

Vt + P{t)xx*h = 0, (6b) 

2 X t + P(t)xl-Mx,t)=0, (6c) 

2a t + \p(t) X „ e -2rf(t)]<r = Q (<r = A,B,C), (6d) 

T t BP T + Z&-rfcCQ m - fiT t CQ + g{t)CQ[C 2 Q 2 + (A + BP) 2 } = 0, (6e) 

-nCQr + ^-^BP m - fiTt (A + BP) + g(t)(A + BP)[C 2 Q 2 + (A + BP) 2 } = 0, (6f) 

where // is a constant, and T](x,t), x(x,i), A(t), B(t), C(t), r(t), P(j],t), Q(t],t) are functions to be deter- 
mined. After some algebra, it follows from Eqs. (|6"aj) - (16dp that we have 

V (x,t) = a(t)x + S(t), x ( a ,t) = - 5 -*L_ s »--^_ a5 + Xb (t), v{x ,t)= Xt + ^- x l, (7a) 
A(t) = a Vl«WI e^ 7(s)ds , B(t) = &A(i), C7(t) = cA(t), (7b) 

where ao, b, c are constants, a(i) (the inverse of the wave width), S(t) (—5(t) /a(t) being the position of its 
center of mass), and xo(^) are &11 f rcc functions of time t. 

To further reduce Eqs. (|6"e|) and (|6fj) to the system of coupled partial differential equations with constant 



coefficients, We require the conditions: r t = ^p-17 2 and g(t) — ^^-GA 2 (t)r] 2 {G = const.), which can generate 
the constraints for the variable r(t) and nonlinearity g(t) 

1 f* 2,^ of ^ Ga(t)(3(t) 



r(t) = - y ^Wd., git) = ^ (8) 



such that Eqs. (IBel) and (|6g reduce to the coupled system of differential equation with constant coefficients 

bP T + cQ m - ficQ + cGQ[c 2 Q 2 + (1 + bP) 2 } = 0, (9a) 
-cQ T + bP m - M(l + W) + G(l + 6P)[c 2 Q 2 + (f + bP) 2 } = 0. (9b) 
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Following the direct approach developed in Refs. [28J129], we can obtain the rational solutions of system (J9J) 
such that the corresponding rational-like solutions (nonautonomous rogons) of Eq. ((2|) can be found in terms of 
similarity transformations ([3]) and In the following, we will exhibit the dynamical behaviors of rational- like 
solutions with many interesting nontrivial features. 

First-order rational-like solution- It follows from system Q that we have the solution P(r), r) = —4/ [bAi (jj, r)] 
and Q{r),r) = —8t/[cAi(t],t)] with A\{t),t) = 1 + 2rj 2 + At 2 for fi = G = 1. Thus based on the similarity 
transformations ([3]) and ([5]), we obtain the first-order rational- like solution (nonautonomous rogon) of Eq. ([2]) 

4 + 8ir(t) 



^(z,«)=a v^Ri)je/° 7(s)ds 
whose intensity can be written as 
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1 + 2[a(t)x + S(t)} 2 + 4r 2 (t) 



e »[x(a!,*)+T(t)] (10) 



|^(x, t )P = agK t )|e 2 /o^^ 2[ ^ + ^ ]2 + 4r2(t) - 3 > 2 + r 2( ^ (ID 
° {l + 2[a(i)a: + <5(t)] 2 +4T 2 (i)} 2 

where x( x >t) an d are given by Eqs. (|7al) and flSJ). 

It is easy to see that the rational-like solution (fT0|) is different from the known rational solution of the NLS 
equation (jTJ) [TSl EH [HI H21 H31 [2H1 since it contains some free functions of time t, which will generate 
abundant structures related to the optical RWs. In particular, when a = 2, ao = /3 = 1, xo — 7 — 0, resulting 
in <7 = 1, Eq. ([2]) reduces to Eq. (flj such that the nonautonomous rogon solution ([TO)) reduces to the known 
rogon solution in Refs. |21[ |2"8] . In what follows, we will choose some free functions of time to exhibit the 
obtained rational-like solution (|T0]) . 

For the fixed parameters ao — 1, Xo(*) = and j(t) = 0.1 tanh(i)sech(t), i) if we choose other free functions 
as the polynomials of time t, i.e. a(t) = 1, Pit) = Q.5t 2 , then Figures la and lb depict the dynamical behavior 
of the rational-like solution (flOl) for different terms 5(t) — t, t 2 , respectively, in which the other coefficients g(t) 
and v(x,t) in Eq. ([2]) are given by 

g( t ) = Ii2 e [ S och( t )-l]/5 ! (12a) 

v(x,t) = Ar^x + r 2 for 5(t)=t, (12b) 
v{x,t) = At- 2 x + A for 5(t) = t 2 ; (12c) 

ii) if we choose other free functions as the periodic functions of time t, i.e. a(t) = dn(t, k), /3{t) — cn(t, k), then 
Figures 2a and 2b display the dynamical behavior of the intensity of the rational-like solution (ITU)) for different 
terms S(t) = sn(t, fc), cn(t, k), respectively, in which the coefficients g(t) and v(x,t) in Eq. ^) are given by 

g(t) = - cn{t, k) dn(i, fc) e t secll (*)- 1 l/ 5 , (13a) 

v(x,t) = k2cn j t,k K 2 + - cn(t,k)[k 2 sd(t,k)x - l} 2 for S(t) = sn(t, *), (13b) 
2dn (t, fc) 2 

u(a;,t)= fc C "^' fc ' ) a: 2 + k \ x + -cn(t, fc)[fc 2 sd(t, k)x + sc(t, fc)] 2 for £(t) = cn(t, fc); (13c) 

2dn (t, k) cn 2 (t, fc) 2 



iii) if we choose other free functions as the periodic functions of time i, i.e. a(t) = cn(i, k), (3(t) — dn(t, k), 
then Figures 3a and 3b exhibit the dynamical behavior of the rational-like solution (fTU|) for different terms 
5(t) = sn(t, k), dn(t, k), respectively, in which the coefficients g(t) and v(x,t) in Eq. © are given by 

g (t) = I cn(i, k) dn(t, fc) e [ soch (*)-i]/5 ; ( 14a ) 

v(x,t) = ^—£pL x 2 + ]-dn(t,k)\sc(t,k)x - l] 2 for S(t) = sn(i,fe), (14b) 
2cn^(i, fe) 2 



u(a;,t)= ^ d "^'^ x 2 + fc C ^ t,k K + ^dn(t,k)\sc(t,k)x + k 2 sd(t,k)} 2 for 5(f) = dn(i, k). (14c) 
2cn^(i, fc) dn (t, k) 2 



It follows from these figures that the rational-like solution (|10l) is different from the known rational solution of 
the NLS equation ([T]) 15, 16, 21, 22, 23, 28, 29 , and may be useful to raise the possibility of relative experiments 
and potential applications for the RW phenomenon. 

Second-order rational-like solution — It follows from Q that we have the second-order rational solution of 
Eq. © in the form P(t),t) = V 2 (r], T)/[bA 2 (r), r)] and Q(v,t) = Q 2 (r),T)/[cA 2 (r),T)}, where 



V 2 (rj(x, t), r(t)) = -\^- Wr 2 10r 4 - \n 2 - 9r 2 + |, 



Q 2 (r,(x,t),T(t)) 



^ + Arfr 2 + 4r 4 - 3r/ 2 + 2r 2 - — 



(15) 



AMx, t), r(t)) = lr, 6 + ir, 4 r 2 + r, 2 r 4 + V + ^ 4 + V - |„V + Ar, 2 + |r 2 + |, 



As a consequence, based on the transformations Q and ([5]), we have the second-order rational-like solution 
(two-rogon solution) of Eq. ([2]) 

V 2 {-n{x,t),r{t))+iQ 2 {-n{x,t),T{t))' 



Mx,t) = a ^(tj\ef°^ ds 
whose intensity is given by 



A 2 (n(x,t),T(t)) 



e i[x(x,t)+T(t)\ (16) 



\M*,t)\ -a Ht)\eJo _____ . 



(17) 



where rj(x,t) = a[t)x + S(t), and x(x, i), r(t) are given by Eqs. (|7aj) and ([8]). The obtained second-order 
rational- like solution (|T6]) is also different from the known rational solutions of the NLS equation (QJ (2TJ [28] , 
since some free functions of time t are involved. Similar to the first-order rational-like solution, the solution 
(fT6|) can also reduce to the known one of Eq. (pj in Refs. [2T| [28]. 

For the fixed parameters ao = 1, Xo(t) = 0, j(t) = 0.1 tanh(i)scch(<), i) if we choose other free functions as 
the polynomials of time t, i.e. a(t) = 1, f3(t) = 0.5t 2 , then Figures 4a and 4b depict the dynamical interaction of 
the two-rogon solution (JTHJ) for different terms 6(t) = O.li, f 2 ; ii) if we choose the free functions as the periodic 
functions of time i, i.e. a(t) — dn(i, k), f3(t) — cn(t, k), then Figures 5a and 5b illustrate the dynamical 
interaction of the two-rogon solution (fTH]) for different terms 6(t) = sn(t, k), cn(t,k); iii) if we take other free 
functions as the periodic functions of time t, i.e. a(t) = cn(i, k), j3(t) — dn(f, k), then Figures 6a and 6b exhibit 
the dynamical interaction of the two-rogon solution (|16j) for different terms 5(t) = sn(f, k), dn(i, k). 
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In conclusion, we have presented the analytical first-order and second-order rational-like solution pairs 
(nonautonomous rogons) of the inhomogeneous nonlinear Schrodinger equation with variable coefficients by 
using the similarity transformation and direct ansatz. By using the direct approach |28) . we can also obtain 
the higher-order rational-like solutions of Eq. ^ which are omitted here. These obtained solutions be used 
to describe the possible formation mechanisms for the optical rogue wave phenomenon in optical fibres, the 
oceanic rogue wave phenomenon in the deep ocean, and the matter rogue wave phenomenon in Bose-Einstein 
condensates. Furthermore, it should be emphasized that we give the explicit solutions for the existence of 
matter rogons in Bose-Einstein condensates. Moreover, the snake propagation traces and the interaction of 
optical nonautonomous rogons are exhibited by choosing some free functions of time. Some shapes of the 
optical nonautonomous rogons and fascinating interactions between two optical nonautonomous rogons are also 
achieved with different functions. These solutions modify the known solutions related to the optical rogons. 
Moreover, this constructive idea can be also extended to other nonlinear systems with variable coefficients to 
generate nonautonomous rogons. This will further excite the study of rogons in the field of nonlinear science. 

Acknowledgements 

This work was supported by the FCT SFRH/BPD/41367/2007 and the NSFC60821002/F02. 

References 

[1] C. Sulem, P. L. Sulem, The Nonlinear Schrodinger Equation: Self-focusing and Wave Collapse, Springer- Verlag, 
New York, 1999. 

[2] Y. S. Kivshar, G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals, Academic, San Diego, 2003. 
[3] A. Hasegawa, M. Matsumoto, Optical Solitons in Fibers, Springer, Berlin, 2003. 

[4] C. J. Pethick, H. Smith, Bose-Einstein Condensation in Dilute Gases, Cambridge University Press, Cambridge, 
2002. 

[5] F. Dalfovo, S. Giorgini, L.P. Pitaevskii, S. Stringari, Rev. Mod. Phys. 71 (1999) 463. 

[6] R. Carretero-Gonzalez, D. J. Frantzeskakis, P. G. Kevrekidis, Nonlinearity, 21 (2008) R139. 

[7] M. Onorato, A. R. Osborne, M. Serio, S. Bertone, Phys. Rev. Lett. 86 (2001) 5831. 

[8] C. Kharif, E. Pelinovsky, Eur. J. Mech. B (Fluids) 22 (2003) 603. 

[9] A. R. Osborne, Nonlinear Ocean Waves, Academic Press, New York, 2009. 

[10] C. Kharif, E. Pelinovsky, A. Slunyaev, Rogue Waves in the Ocean, Observation, Theories and Modeling, Springer, 
New York, 2009. 

[11] L. Draper, Mar. Obs. 35 (1965) 193. 

[12] P. Miiller, Ch. Garrett, A. Osborne, Oceanography 18 (2005) 66. 

[13] D. R. Solli, C. Ropers, P. Koonath, B. Jalali, Nature 450 (2007) 1054. 

[14] D.-I. Yeom, B. Eggleton, Nature 450 (2007) 953. 

[15] D. H. Peregrine, J. Aust. Math. Soc. Ser. B, Appl. Math. 25 (1983) 16. 

[16] Y. C. Ma, Stud. Appl. Math. 60 (1979) 43. 



G 



[17] N. Akhmediev, V. I. Korneev, Theor. Math. Phys. 69 (1986) 1089. 

[18] N. Akhmediev, V. M. Eleonskii, N. E. Kulagin, Theor. Math. Phys. 72 (1987) 809. 

[19] K. B. Dysthe, K. Trulsen, Phys. Scr. T82 (1999) 48. 

[20] V. V. Voronovich, V. I. Shrira, G. Thomas, J. Fluid Mech. 604 (2008) 263. 
[21] N. Akhmediev, J. M. Soto-Crespo, A. Ankiewicz, Phys. Lett. A 373 (2009) 2137. 
[22] N. Akhmediev, J. M. Soto-Crespo, A. Ankiewicz, Phys. Rev. A 80 (2009) 043818. 
[23] A. Ankiewicz, N. Devine, N. Akhmediev, Phys. Lett. A 373 (2009) 3997. 
[24] M. Onorato, A. R. Osborne, M. Serio, Phys. Rev. Lett. 96 (2006) 014503. 

[25] P. K. Shukla, I. Kourakis, B. Eliasson, M. Markhmd, L. Stenflo, Phys. Rev. Lett. 97 (2006) 094501. 
[26] A. Gronlund, B. Eliasson, M. Markhmd, Eur. Phys. Lett. 86 (2009) 24001. 
[27] V. P. Ruban, Phys. Rev. Lett. 99 (2007) 044502. 

[28] N. Akhmediev, J. M. Soto-Crespo, A. Ankiewicz, Phys. Rev. E 80 (2009) 026601. 
[29] N. Akhmediev, A. Ankiewicz, M. Taki, Phys. Lett. A 373 (2009) 675. 
[30] M. Hopkin, Nature 430 (2004) 492. 

[31] D. R. Solli, C. Ropers, B. Jalali, Phys. Rev. Lett. 101 (2008) 233902. 
[32] K. Hammani, C. Finot, J. M. Dudley, G. Millotl, Opt. Exp. 16 (2008) 16467. 
[33] J. Kasparian, P. Bejot, J.-P. Wolf, J. M. Dudley, Opt. Exp. 17 (2009) 12070. 
[34] Yu. V. Bludov, V. V. Konotop, N. Akhmediev, Phys. Rev. A 80 (2009) 033610. 
[35] B. A. Malomed, J. Opt. B 7 (2005) R53. 

[36] S. Ponomarenko, G. P. Agrawal, Phys. Rev. Lett. 97 (2006) 13901. 

[37] V. N. Serkin, A. Hasegawa, T. L. Belyaeva, Phys. Rev. Lett. 98 (2007) 074102. 

[38] M. Centurion, M. A. Porter, P. G. Kevrekidis, D. Psaltis, Phys. Rev. Lett. 97 (2006) 033903. 

[39] V. M. Perez-Garcia, P. J. Torres, V. V. Konotop, Physica D 221 (2006) 31. 

[40] J. Belmonte-Beitia, V. M. Perez-Garcia, V. Vekslerchik, V. V. Konotop, Phys. Rev. Lett. 98 (2008) 064102. 
[41] Z. Y. Yan, Phys. Lett. A 361 (2007) 223. 
[42] Z. Y. Yan, Phys. Scr. 75 (2007) 320. 

[43] Z. Y. Yan, Constructive Theory and Applications of Complex Nonlinear Waves, Science Press, Beijing, 2007. 
[44] Z. Y. Yan, V. V. Konotop, Phys. Rev. E 80 (2009) 036607. 
[45] N. J. Zabusky, M. D. Kruskal, Phys. Rev. Lett. 15 (1965) 240. 



7 



List of the Figure Captions 



Figure 1. Wave propagations (left column ) and contour plots (right column) for the intensity |V'i| 2 (fTTj) of the first- 
order rational-like solution (jlOl for a = o.q = 1.0, f3 — 0.5£ 2 , 7(<) = 0.1 tanh(i)sech(i). (a)-(b) 8(t) = t; 
(c)-(d) 5{t)=t 2 . 

Figure 2. Wave propagations (left column) and contour plots (right column) for the intensity \ipi\ 2 (fTTj) of the first- 
order rational-like solution (fT0|) for ao = 1-0, j(t) — 0.1 tanh(i)sech(i), k = 0.6, a = dn(i, fc), (3(t) = 
cn(t,fe): (a)-(b) 5{t) =sn(t,k); (c)-(d) 5{t) =cn(t,k). 

Figure 3. Wave propagations (left column) and contour plots (right column) for the intensity \ipi\ 2 (fTTj) of the first- 
order rational-like solution (fT0|) for «o = 1.0, j(t) = 0.1tanh(t)sech(t), k — 0.6, a = cn(t,k), j3(t) = 
dn(t,fc): (a)-(b) 5{t) =sn(i,*); (c)-(d) 6(t) =dn(t,A). 

Figure 4. Wave propagations (left column) and contour plots (right column) for the intensity \4>2\ 2 (|T7)) of the second- 
order rational-like solution (IT^l) for a = 1-0, 7(i) = 0.1tanh(t)sech(i), k = 0.6, a = dn(t,k), (3(t) = 
cn(t,k): (a)-(b) 5(t) = O.lt; (c)-(d) £(t) = t 2 . 

Figure 5. Wave propagations (left column) and contour plots (right column) for the intensity V^l 2 (fl7|) of the second- 
order rational-like solution (fl~6"l) for ao = 1.0, 7(f) = 0.1tanh(t)sech(i), fc = 0.6, a = dn(i, fe), (3(t) = 
cn(t,k): (a)-(b) 5(t) =sn(*,fc); (c)-(d) S(t) = cn(t,fc). 

Figure 6. Wave propagations (left column) and contour plots (right column) for the intensity V^l 2 (|17[) of the second- 
order rational-like solution (fT6|) for ao = 1.0, j(t) = 0.1 tanh(i)sech(i), fc = 0.6, a = cn(t,k), f3(t) = 
dn(t,fc): (a)-(b) £(t) = sn(t,fc); (c)-(d) *(f) =dn(t,fc). 




Figure 1: (color online). Wave propagations (left column) and contour plots (right column) for the intensity l^il 2 

of the first-order rational-like solution (JTUJ for a = ctQ = 1.0, = 0.5i 2 , j(t) = 0.1tanh(t)sech(t). (a)-(b) 5{t) = t; 




Figure 2: (color online). Wave propagations (left column) and contour plots (right column) for the intensity l^il 2 (1 1 1 1> of 
the first-order rational-like solution (|10[1 for ao — 1.0, fit) = 0.1 tanh(t)sech(f), k = 0.6, a(t) = dn(t, k), fi{t) = cn(t, k): 
(a)-(b) S(t) = sn(t, fe); (c)-(d) S(t) = cn(t, fc). 
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Figure 3: (color online). Wave propagations (left column) and contour plots (right column) for the intensity \ipi\ 2 (1111) of 
the first-order rational-like solution (|10[) for a?o = 1-0, y(i) = 0.1 tanh(t)sech(£), k = 0.6, a(t) = cn(t,k), /3(t) = dn(t,k): 
(a)-(b) 5(t) = sn(t, k); (c)-(d) S(t) = dn(t, fe). 




Figure 4: (color online). Wave propagations (left column) and contour plots (right column) for the intensity \ip2\ 2 (1 1 7ft of 
the second-order rational-like solution 1)161) for a = oto = 1.0, 7(t) = 0.1tanh(t)sech(t), /3(t) = 0.5t 2 : (a)-(b) S(t) = O.lt; 
(c)-(d)5(t)=t 2 . 
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Figure 5: (color online). Wave propagations (left column) and contour plots (right column) for the intensity \ip2\ 2 (|17[) 
of the second-order rational-like solution (|16p for «o — 1-0, j(t) = 0.1 tanh(t)sech(£), k — 0.6, a(t) — dn(t,k), /3(t) = 
cn(t, fe): (a)-(b) 5(t) = sn(t, k); (c)-(d) 5(t) = cn(t, k). 




Figure 6: (color online). Wave propagations (left column) and contour plots (right column) for the intensity \ip2\ 2 (|17[) 
of the second-order rational-like solution (|16|l for Qo = 1.0, ^(t) — 0.1 tanh(f)sech(t), k = 0.6, a(t) = cn(t,k), f3(t) = 
dn(t, fc): (a)-(b) <5(t) = sn(i, fc); (c)-(d) 5(t) = dn(t, fc). 
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